Colloidal transport through optical tweezer arrays 
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Viscously damped particles driven past an evenly spaced array of potential energy wells or barriers 
may become kinetically locked in to the array, or else may escape from the array. The transition be- 
tween locked-in and free-running states has been predicted to depend sensitively on the ratio between 
the particles' size and the separation between wells. This prediction is confirmed by measurements 
on monodisperse colloidal spheres driven through arrays of holographic optical traps. 



Particles rolling down a sinusoidally modulated slope 
constitute an archetype for problems as diverse as the 
electrodiffusion of atoms on crystals, the transport char- 
acteristics of Josephson junctions and the entire field of 
chemical kinetics. The one-dimensional tilted washboard 
problem has been studied exhaustively. Considerably less 
is known for related problems in higher dimensions. Re- 
cently, attention has become focused on the transport 
of viscously damped colloidal particles flowing through 
two-dimensional potential energy landscapes. This sys- 
tem can be realized in practice by passing fluid-borne ob- 
jects through microfabricated arrays of posts QII2>OJ0> 
over arrays of electrodes JSH and through periodically 
structured light fields 0, U Hi- 
Experimental realizations of two-dimensionally modu- 
lated transport are interesting both because they provide 




FIG. 1: Colloidal transport through an array of holographic 
optical traps, (a) Holographic trapping system, (b) Phase 
hologram <p(p) encoding at 20x2 array of optical traps, shown 
in (c). Scale bar indicates 5 fim. 



insights into the underlying fundamental problem, and 
also because they constitute an entirely new category of 
sorting techniques. Different types of objects, it turns 
out, can follow dramatically different paths through the 
same physical landscape. Sorting fluid-borne objects by 
size, shape and composition have been demonstrated in 
this way [1,13 ■ Preliminary theoretical studies 0,^1 SU S~ 
gest that periodic landscapes can act as extraordinarily 
selective sieves, for example sorting spheres by size with 
exponential resolution. This article presents experimen- 
tal confirmation of some of these theories' predictions. 

Our system, shown schematically in Fig. ^ tracks 
the motions of monodisperse colloidal spheres as they 
are driven back and forth over static potential energy 
landscapes created with arrays of holographic optical 
traps. The samples consist of colloidal silica spheres 
a = 0.75 ±0.075 ^m in radius (Bangs Laboratories 5303) 
These spheres were dispersed in deionized water and her- 
metically sealed in a slit pore formed by bonding the 
edges of a #1 cover slip to the face of a glass micro- 
scope slide. The glass surfaces were treated by oxygen 
plasma etching before assembly to increase their surface 
charge and thereby prevent particle deposition. The sam- 
ple was rigidly mounted on a Prior ProScan II translation 
stage integrated into a Nikon TE2000U optical micro- 
scope, where it was allowed to equilibrate to room tem- 
perature, T = 296 ± 2 K. 

Previous experimental studies of transport through 
static light fields have driven the spheres electrokineti- 
cally, optophoretically or hydrodynamically 0, @, 01 ■ 
Others have swept the l ight field through stationary sam- 
ples [H [H Q |H H[]j. We instead used the motor- 
ized stage to translate the entire sample past stationary 
patterns of optical traps. All particles consequently trav- 
eled past the traps at the same velocity, v, without com- 
plications due to nonuniform flow profiles |9j and without 
time-dependent ratchet phenomena 0|. Roughly 5,000 
particles were repeatedly passed back and forth over the 
same field of view at constant speed and a variety of an- 
gles to build up a statistically well sampled set of data 
for each set of conditions. Repeatedly revisiting the same 
part of the sample cell with the same particles minimized 
effects due to nonuniform sample thickness and variabil- 
ity in the spheres' properties. 

Images of the moving particles were recorded as an un- 
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FIG. 2: Selected trajectories of a — 0.75 fim diameter silica 
spheres flowing at v = 10 /jm/s through an array of optical 
traps with b — 2.025 /im at P = 0.85 W. The shaded boxes 
indicate the extent of the optical trap array, (a) Aligned: 
9 = 0°. (b) Kinetically locked in at 9 = 10°. (c) Marginally 
locked in at 9 m — 12°. For 9 > 8 m , some particles escape the 
array. 

compressed digital video stream using an NEC TI-324A 
video camera and a Pioneer 520HS digital video recorder. 
The combination of a 100 x oil immersion objective lens 
(Nikon Plan Apo, NA 1.4) and a 0.63x video eyepiece 
provides a field of view of 86 x 63 ^tm 2 and a magnifi- 
cation of 0.135 /xm per pixel. Spheres' images were sub- 
sequently digitized into trajectories with 20 nm spatial 
resolution at 1/30 s intervals using standard methods of 
digital video microscopy |l8j |. Typical measured trajec- 
tories are plotted in Fig. [21 

The moving spheres encountered point-like optical 
tweezers [l9j projected in a 20 x 2 gri d with the holo- 
graphic optical trapping technique |20l l2ll |22| . All of 
the traps were created from a single beam of laser light 
at wavelength A = 532 nm and power P = 0.85 W pro- 
vided by a frequency-doubled diode pumped solid state 
laser (Coherent Verdi). This beam was diffracted into an 
array of trap-forming beams by a phase-only computer- 
generated hologram, </?(p), an example of which ap- 
pears in Fig. ^b). The diffracted beams then were 
relayed to the objective lens, which focused them into 
traps. Imprinting the hologram on the beam's wave- 
fronts with a computer-addressed spatial light modula- 
tor (SLM, Hamamatsu X7269 PPM) allows for sequences 
of trapping patterns to be projected with different lat- 
tice constants and orientations. Laser light is reflected 
into the objective lens with a tuned dichroic mirror 
(Chroma Technologies) with a reflectivity of 99.5 per- 



cent at A = 532 nm. This mirror transmits light at other 
wavelengths, which therefore can be used to create im- 
ages of the spheres. 

The image of the focused traps in Fig. ^c) was ob- 
tained by placing a front surface mirror in the lens' focal 
plane. Enough of the reflected laser light passes through 
the dichroic mirror to create a clear low-noise image of 
the trapping pattern. Images such as this were used to 
adaptively improve the traps' uniformity |22j |. After iter- 
ative improvement, the traps' intensities typically varied 
by less than 15 percent from the mean. 

Each colloidal particle experiences an optical trap as 
radially symmetric potential energy well whose depth, 
Vo(a), and width, a(a), both depend on the particle's 
shape, size and composition ,10]. The trap's depth also 
is proportional to the total laser power, P. If a particle 
is deflected enough by its encounter with one trap to fall 
into another trap's domain of influence, it can become ki- 
netically locked in to a commensurate trajectory through 
the array of traps 0,&EH13- If> on the other hand, the 
driving force is too strong, or the required deflection an- 
gle too steep, the particle escapes from the traps and runs 
freely downstream. In our experiment, the driving force, 
F(j(a) = 7(a) v, is the hydrodynamic drag on a sphere of 
radius a driven through a quiescent fluid at velocity v. 
The viscous drag coefficient, 7(a), is proportional to the 
fluid's viscosity and accounts for hydrodynamic coupling 
to the bounding walls pi I2I I2I w\ . 

The data in Fig. were obtained at v = 10 /im/s so 
that each sphere crossed the entire field of view within 
t = 9 s. By contrast, the spheres' measured 0, ^ 
self-diffusion coefficient, D = 0.10 ± 0.01 /im 2 /s, cor- 
responds to a thermally driven displacement of just 
\/2Dt = 0.6 fim in the same period. The associated 
viscous drag coefficient, 7 = ksT/D — 40 fN s//xm sug- 
gests that the spheres were driven past the traps with a 
maximum force of roughly Fq = 0.4 pN. 

The distinction between locked-in and free trajectories 
becomes clear when the driving force Fq is inclined with 
respect to the array, as shown in Fig.^b). In this case, 
the locked-in trajectories are deflected by angle 8 with 
respect to Fq. This deflection is the basis for continuous 
sorting techniques in which different fractions of a mixed 
sample are deflected to different angles by the same op- 
tical intensity field 0,HI3|. 

All trajectories become locked in when the driving 
force is aligned with the traps, 8 = 0°; they all escape 
for 8 = 90°. The maximum angle, 8 m , to which a spheri- 
cal object can be deflected by a periodic optical intensity 
field before it escapes has been predicted to depend ex- 
ceptionally strongly on the object's radius 0, 0, l28| . 
This is a purely kinematic effect, and not the result of 
thermal activation over potential energy barriers, as has 
been suggested @- 

Modeling the optical traps as Gaussian potential en- 
ergy wells separated by distance b yields |9lll0| 

sinfl m = \S(a)\ exp (- \ - Y (1) 
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with the prefactor, 



2 V (a) 



a(a)F (a) ' 



(2) 



reflecting a balance between trapping and driving forces 
at the point of escape. The traps' effective width is given 
by 



A 2 

a 2 (a) = a 2 + 
w An 2 



(3) 



where n is the fluid's index of refraction. For the par- 
ticles in our study, dispersed in water with n = 1.33, 
a = 0.85 /xm. Similar results are obtained for more gen- 
eral periodic landscapes [Hj, with Eq. JJJ representing 
the leading term in a Fourier expansion of the trapping 
potential. Equation Q also applies for arrays of barriers 
rather than wells jlOj . 

According to Eq. the marginally locked-in angle, 
9 mi depends very strongly on the ratio b/a(a), and thus 
on particle size for a > A. This is the basis for the asser- 
tion |9j that sorting by transport through a periodic land- 
scape can offer exponential size selectivity. This predic- 
tion, however, results from limiting arguments that have 
not been tested experimentally. More rigorous results 
are available only in the limit that Vq <C F a f=a fcgT, in 
which case thermal forces cannot be neglected and the 
escape transition is expected to be less dramatic 29] . 

We explicitly tested the predicted dependence on trap 
separation, b, and laser power by driving monodisperse 
particles past adaptively optimized trap arrays at fixed 
speed, v, over a range of angles, gauging the marginal 
angle by the suddenly increasing proportion of escaping 
trajectories. As 6 approaches 8 m , particles' trajectories 
become increasingly sensitive to variations in the traps' 
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FIG. 3: Dependence of the marginally locked in angle as a 
function of lattice constant. Inset: The data replotted to 
emphasize the comparison to Eq. Q. 
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FIG. 4: Dependence of the marginally locked in angle on 
total projected laser power for fixed values of lattice constant, 
b — 1.35 fj,m, and speed, v — 46 /im/s. The square point 
is replotted from Fig. |H] Inset: Probability density p(r) to 
find a particle at r, normalized by the bulk value, po, at 6 = 
45° . The downstream shadow in p(r) demonstrates effectively 
perfect large-angle deflection for P = 1.5 W. Average over 
500 trajectories. 



intensities, to thermal fluctuations, and to small differ- 
ences in individual particles' radii. A typical escape tran- 
sition is apparent in the selected trajectories plotted in 
Fig-EJ c ) . It is reasonable to expect that the first observed 
escape events would involve the smallest particles inter- 
acting with the weakest traps. Analyzing particle tracks, 
however, did not provide sufficient size resolution to test 
this directly. The marginally locked-in angle for a par- 
ticular choice of b was determined by analyzing roughly 
2000 such trajectories for each value of 6 ranging from 0° 
to 30° in 2° increments. 

Results for 9 m (b), obtained at fixed speed, v, and laser 
power, P, are plotted in Fig. |3J The dashed curve is 
a one-parameter fit to Eq. Q, where only the overall 
scale, S(a), is treated as a free parameter. The fit value, 
S(a) = 0.40±0.01 corresponds to a well depth of roughly 
100 ksT, which is reasonable for traps powered with 
2 mW of light. The experimental results' excellent agree- 
ment with the prediction suggests that Eq. JIJ quantita- 
tively describes colloidal transport through an array of 
optical traps. 

The marginally locked-in angle can be tuned for maxi- 
mum sensitivity by adjusting b. Its magnitude can be set 
independently through the prefactor, S(a). In particular, 
S(a) should scale linearly with laser power, P, through 
its dependence on Vq. Thus 8 m for a given size of particle 
can be adjusted with laser power, for fixed v, up to the 
point that particles begin to get stuck in the traps. The 
data in Fig. ^demonstrate that deflection angles as large 
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as 6 m = 45° can be attained in this way. They therefore 
contradict the assertion ||| that arrays of discrete traps 
are incapable of deflecting trajectories over large angles. 

The data in Fig. 0] were obtained with a single row of 
20 traps at fixed inter-trap separation of b — 1.35 /jm, 
and with particles driven at v = 46 /im/s. Using a single 
line of traps increases the possibility that particles might 
leak through the array, but doubles the accessible range 
of laser powers. Data from Fig. [3] at the same particle 
speed and lattice constant fall on the same curve once 
the laser power is rescaled. The observed linear depen- 
dence of sm# m on laser power, P, confirms the predicted 
dependence on the prefactor, S(a) in Eq. QJ. 

The plot inset into Fig. 01 provides an overview of the 
data obtained at 9 m = 45° and P = 1.5 W. It shows 
the probability distribution p(r) for finding a particle 
within 200 nm of r, integrated over 500 trajectories, 
comprising roughly 30,000 separate particle images. The 
probability distribution is normalized by the probability, 
poj °f finding a particle in the undeflected stream over 
the same period. Darker regions indicate a higher-than- 
average probability density and lighter regions indicate 
lower-than-average probabilities. Details in the distribu- 
tion are emphasized with a nonlinear color table, which 



also is inset. The results show that particles are strongly 
concentrated in the traps themselves and work their way 
up the array until they escape at its end. The nearly per- 
fect deflection of the incident stream of particles leaves a 
sharply defined shadow in the probability density down- 
stream of the traps. A total of 15 of the 500 particles 
escaped the array near its end, most likely representing 
the smallest end of the particle size distribution. 

Confirming the marginally locked-in angle's depen- 
dence on trap separation and laser power supports the 
assumptions made in deriving Eq. l[T]l. particularly be- 
cause the predicted form for <j(a) agrees quantitatively 
with experimental results. This provides additional, al- 
beit indirect, support for the prediction that athermal 
sieving by periodically modulated landscapes can sort 
objects with exponential size selectivity. It leaves open 
questions regarding the nature of transitions among dif- 
ferent commensurate locked-in states in more extensive 
two-dimensional lattices. It also does not address the na- 
ture of colloidal transport through aperiodic landscapes, 
such as quasiperiodic arrays of optical traps. Experi- 
ments to address these questions are in progress. 

This work was supported by the National Science 
Foundation through Grant Number DMR-0451589. 
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